Instability of charged and rotating naked singularities 
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We provide evidence that "super-extremal" black hole space-times (either with charge larger 
than mass or angular momentum larger than mass), which contain naked singularities, are unstable 
under linearized perturbations. This is given by an infinite family of exact unstable solutions in 
the charged non rotating case, and by a set of (unstable) numerical solutions in the rotating case. 
These results may be relevant to the expectation that these space-times cannot be the endpoint of 
physical gravitational collapse. 



It is well known that the equations of general relativity admit exact solutions that contain singularities. In some 
• cases, as in black hole space-times, the singularities are "hidden" behind horizons, that is, regions that cannot 
communicate causally with the rest of space-time. However, the same solutions that describe black holes can, for 
' certain choices of their parameters, describe "naked" singularities. This is the case of a charged (Reissner-Nordstrom) 
^ \ space-time with \Q\ > M or a rotating (Kerr) solution with angular momentum a > M . These represent problematic 
D ■ space-times, since the singularities can communicate causally with the exterior. There has always been the expectation 
\ that these solutions are not the endpoint of physically relevant collapsing matter. This whole issue frames itself in the 

■ context of the "cosmic censorship" hypothesis of Penrose that loosely stated claims no naked singularities are formed 
I by gravitational collapse. Several arguments have been put forward that show that one indeed cannot drop charges 

■ into a charged black hole to turn it into "super-extremality" or spin up a Kerr black hole into a similar regime [ij. 
^ ' In this paper we would like to argue that "super-extremal" space-times are unstable under linearized perturbations. 

This would further strengthen the idea that they cannot be the endpoint of gravitational collapse. We start our 
discussion by considering a spherically symmetric, charged space-time described by the Reissner-Nordstrom solution 
. of the Einstein equations, 

o ■ 

\C>\ ds^ = -fdf + {l/f)dr^ + r'^dV? (1) 



O 



where f — {r^ — 2rM + Q'^)/r'^, and dil^ is the standard line element on the unit sphere. The real (positive or 
qh, negative) constants M and Q can be respectively identified with the Newtonian mass and the total (electric) charge 
of the source. We consider only the case M > in what follows. The case M < 0, Q = has been previusly analyzed 
in and 

The range of r in ([1]) is limited by the singularities of the metric coefficients. These occur for ?' and r± — 
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M ± M"^ — Q^. For M > \Q\ the singularity at r = r+ corresponds to a regular horizon, and the solution may be 
extended to r < r+ . The resulting space-time is known as a "charged black hole" . The horizon hides a curvature 
singularity for r ~ 0, rendering the portion of space time outside the horizon globally hyperbolic. However, if \Q\ > M, 
there are no horizons, and the space-time contains a "naked" singularity for r = 0. 

The question of the (linear) stability of both charged and uncharged black hole space-times has been analyzed 
by many authors, starting with the papers by Regge and Wheeler, and Zcrilli fi| (see [HI for a recent formulation 
in arbitrary dimensions). To simplify the discussion we shall restrict to polar (scalar) metric perturbations in the 
Regge- Wheeler gauge [j, @] ■ The metric and electromagnetic scalar perturbations for the angular modes with fixed 
£, {£ = 2,3, ...), can then be encoded in two functions, ^i{r,t),i — 1,2, that satisfy the equations 

'""^ ^/3#+AV.^ + ./.V.^-A<I>. (2) 



dt"^ dr*^ \ dr* 
where k = {i - !)£{£ + 1)(£ + 2), (3i = 3(M -I- (-l)*Ai), and 



with fj, = y/M'^ + (£ -!){£ + 2)/9, and r* the "tortoise" coordinate defined by dr*/dr = 1//, 

/r2-2Mr + g2\ 2KP - ( r~M \ 

r* = r + M\n — — + , arctan ^ + c. (4) 
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These have been analyzed in the case > , where there is a regular horizon, hiding the curvature singularity 
at r = 0. On the other hand, we are here interested in the case > M^, with M > 0. In this case, the range of 
r is the full interval (0, +oo), so it is crucial to establish what boundary conditions are acceptable for r ^ 0. If we 
assume that the solutions are regular near r ~ 0, and expand in power series in r, with time dependent coefficients, 
replacing these expansions in ([2]), we find. 



{t, r) ^ b, (t) r + c^ it) + 



[2{i + 2){l - 1)Q^ + 3Af (Af - Mc, jt) 

8Q4 + Q2 



0{ 



(5) 



where 61.2, ci_2 are, in principle, arbitrary functions of t. Notice that both $i contain two arbitrary functions of 
t and, therefore, these expansions represent the general behavior of these functions near r — 0. We may replace 
the expansions in the expressions for the perturbations of the metric in the Regge- Wheeler gauge near r — 0. It 
is then straightforward to obtain, for instance, the Kretschmann invariant, /C = RabcdR"^'"^'^ up to first order in the 
perturbations. We recall that for the background metric we have K. = 56Q*/r^ + ©(r"^). The explicit expression 
including first order terms is rather long to give explicitly here. After replacing the metric coefficients in /C and some 
simplification we find that, in general, the perturbations introduce singular contributions that diverge faster than the 
background, unless the following conditions are imposed 



Ci(t) 



3 (M + fi) 



blit) C2{t) = 



3 (M - ^) 
4Q2 



hit) 



(6) 



in which case the Kretschmann invariant has the behavior given above. We shall take ([6|) as the appropriate boundary 
condition for the perturbations in what follows. We therefore exclude perturbations that might either modify the 
nature of the singularity, or take us beyond the domain of linear perturbations. We remark that these are the most 
restrictive boundary conditions that can be imposed at the singularity without making the solution trivial. We will 
see, nevertheless, that even under this restrictions the perturbation equations have unstable solutions. 

The usual procedure for analyzing stability is to consider solutions of Q of the form $i(t,r) = eKp{—iLut)(j)i{r) and 
look for pure imaginary values of lo. Replacement of this Ansatz in ^ leads to second order ordinary differential 
equations for 0^ . For the negative mass Schwarzschild singularity, analytic expressions for unstable modes of arbitrary 

satisfying appropriate boundary conditions, were obtained in [Sj]- R was then noticed by Cardoso and Cavaglia 
Q that these unstable modes agree with the algebraically special (AS) solutions in 0], which, although irrelevant as 
perturbations in the black hole (Af > 0) regime -due to their behavior at the horizon- satisfy appropriate boundary 
conditions in the nakedly singular case. The AS modes for the Reissner-Nordstrom spacetime have pure imaginary 
frequency lu — —ik with k — — k/(2/3i) (so that k > 0). They can be constructed following Q, the result being 



$2=0, 



$1 = 



Cixir) + C2x{r) 



du 



kt 



where 



Xir) 



-kr* 



{i-l){e + 2)r + Pi 



with k = 



{£ -!)£{£+ l)(e + 2) 



(v/9Af2 + 4Q2(^_l)(^ + 2) - 3 A/) 



(7) 



(8) 



We can readily check that for all > M^, and i, unstable solutions satisfying appropriate boundary conditions are 
obtained by simply setting C2 = in ([7]). Nevertheless, since /3i < 0, the resulting solutions for $1 are singular for 
r = —l3i/[{£ — + 2)]. We notice, however, that replacing these expressions into those for the metric perturbations 
we find that this "kinematic singularity" is absent both in the electromagnetic field and the perturbations of the metric 
coefficients, and, therefore, for finite t these solutions correspond to the evolution of regular perturbations that can 
be made initially arbitrarily small as compared to the background (see Q for a discussion of an analogous situation 
in the Q = 0, M < case). This is confirmed by an explicit computation of /C to first order in the perturbations. The 
result (exact in r) is, 



_ 8(6r2Af 2 _ 12rAf g2 + 7q4) ce {£ +l){£ + 2) (£ - 1) {Mr - Q^) e^'e-'^^y^™ 

where the first term on the R.H.S of (O corresponds to the background, and C is an arbitrary constant. We therefore 
see that for all > M'^, there exist divergent perturbations of the Reissner-Nordstrom spacetime for all £ > 2. 
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Notice that if < M^, the exp(— fcr*) factor in ([8]) gives a divergent behavior at the horizon for any choice of 



Ci, since r* ~ Mln(r — r"*") if r 



,r > 



and fc > 0, and, therefore, ([8]) cannot be considered an acceptable 



perturbation of the Reissner-Nordstrom spacetime. 

The issue of dynamics in non-globally-hyperboHc static space-times is a subtle one, some aspects of which were 
analyzed in a sequel of papers by Wald |8|] and Wald and Ishibashi 0, [l^]- In d, Q the Klein Gordon equation 
(with positive mass) is cast into the form A similar form for gravitational perturbations on AdS backgrounds 
of arbitrary dimensions is given in [l^]- In all cases, the fact that the analogous of the operator Ai in eq. ([2]) is 
positive definite on an appropriate function space, allows to globally define the dynamics of the field (see, e.g., eq.(4) 
in A straightforward consequence of the existence of the exponentially growing modes Q-jH) is, however, that 

Ai in equation ([2|) is not positive definite on the function space of physically relevant perturbations, i.e., those that 
do not increase the degree of divergence of the scalar invariants as the singularity is approached. Thus, gravitational 
perturbations around super extremal charged black holes lie outside the scope of d, U [l^ . A similar situation was 
found in the non-globally- hyperbolic negative mass Schwarzschild solution, and is discussed in detail in [3|. 

We now turn our attention to the uncharged, rotating Kerr space-times with a > M. The linearized perturbations 
were first studied by Teukolsky [nj, who showed that they could be captured in a master equation in terms of a lin- 
earized tetrad component ip of the Weyl tensor. The equation can be separated by assuming ip = F(r, 9, tp) exp(— ztji) 
and F(r, 9, ip) — S^{6) exp(zr7i(^)i?„,f^m(r), which leads to a coupled system for S and R, 



1 d f . dS 
sin 9 d9 V'' d9 



c?u? cos^ ( 



laujs cos 9 



2ms cos ( 



^d'^R ^, ^^,dR ( - 2is(r - M)K 
A^ + 2(. + l)(.-M)-+ i >- 



s^cot^9 + E-s^] 5 = (10) 
Airujs - X)-R = (11) 



where s = ±2, A = - 2Mr 



and K = (r 



a 

-.2, 



iui — am. 
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The eigenvalues E are determined by regularity 
conditions on S{9) for 9 = 0,7t, and \ = E — 2amuj + a^u)'^ — s{s + 1). 

For non extremal Kerr black holes (a^ < M^), where a horizon is present, the region outside the horizon is globally 
hyperbolic, and one can show stability of the perturbations under appropriate boundary conditions at the horizon 
and for r — > oo. However, for > M^, the space-time is not globally hyperbolic, and even though (|10p still applies, 
the question of the stability is certainly more subtle. As a preliminary Ansatz, we will assume that even in this case 
any acceptable unstable perturbation (solution of PH)) and (fTTI) ') must fulfill at least the requirements that i) it can 
be made arbitrarily small at some chosen time, and ii) that it grows exponentially in time. We think that prior to 
the present analysis, it was not known if even these modest requirements could in fact be satisfied. 

As shown in [6[, AS modes oc exp{^iujt + im(j>) exist for the Teukolsky equation if the Starobinsky constant vanishes. 
This condition and the regularity condition S{9 = 0, tt) = required in pop . impose two constraints on A and co that 
may have complex lu solutions These solutions are of the form R(r) = (^ -I- Br + Cr^ + Dr^) exp(ia;r*) where 
r* ^ r as \r\ — > oo, and, since the domain of interest when > is — oo < r < oo, AS modes always diverge in one 
of the asymptotic regions. Therefore, we do not consider AS modes as relevant perturbations of the super-extremal 
Kerr solution because they do not satisfy our Ansatz. 

In general, solutions of the Teukolsky system can be carried out numerically, as first discussed by Press and 
Teukolsky [12]. In our case we started by assuming that lo = ik with /c > to seek for solutions that grow with time. 
For simplicity we set m = 0, and chose E to be real, so the equation for R had real coefficients. The equation for S 
is complex. We integrated (jlip using a power series in x = cos 9 around x = 0, up to the highest order allowed by 
the particular computer implementation that we used, which in our case was x^'', but we checked that essentially the 
same results were already obtained using expansions up to order 20 or higher. Regularity of S for 9 — Q^-k, implies 
S{x = ±1) = 0. Imposing this condition yields a lengthy yet polynomial relationship between E and fc. We chose the 
lowest real value of E given fc. With this value of E we solved numerically equation pT|) using a shooting method. 
To set up the shooting method we worked out asymptotic approximations to the solution for large values of |r| (in 
the super-extremal case Kerr can be extended through the "ring" singularity to negative values of r) . Indicating with 
± the cases r — > ±oo, they are of the form. 



Rir) 



-k\r\-(c^±kM) In(r^) 



(12) 



where c~ = 1/2, = s + 1/2, and the constants are adjusted so that pT|) is satisfied to the given order. These 
expansions were used to generate numerical initial data for the shooting algorithm, typically around r = ±10 (M = 1). 
For a we took values in the range 1.1 — 1.4. The value of k was then varied until the algorithm yielded a finite solution. 
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cos(e) 



FIG. 1: The left panel shows a typical solution for the radial part of the perturbation R{r) in the Kerr case, obtained for 
s = —2, M — 1, and a = 1.4. The amplitude is arbitrary since this is a linear problem, and has been normalized so that 
R{r) = 1 at its maximum. The right panel shows the corresponding real and imaginary parts of the angular function S{6). 
The integration gave k = 7.07..., and E = 17.37... 



A typical solution for s = —2 is shown in figure 1. Similar results, for the same values of fc, were obtained for s = +2, 
(i.e., for "^o) as expected from the non vanishing of the Starobinsky constant [5]. 

The shooting algorithm was implemented in Maple using the built in Runge-Kutta integrator. We checked that the 
solution was insensitive to the error tolerance of the integrator and to the departure point of the shooting. The solution 
has an exponential behavior that limits in practice how far in the range of r can we start the shooting procedure 
before running into machine precision problems. We verified the solution by shooting both from large negative and 
positive values of r. The conclusion of the numerical analysis is that solutions with A: > exist, at least for the values 
of the other parameters chosen, and seem to be numerically robust. Further studies are needed to confirm the ranges 
of values of the various parameters for which unstable solutions exist. 

It should be stressed that any solution i/j of the s = — 2 Teukolsky equation is a "Debye potential" from which a 
metric perturbation hab that solves the linearized Einstein equations around the Kerr background can be constructed. 
hab is obtained by applying a second order linear differential operator to ip (see [isj). Using the solution above as 
a Debye potential gives a linearized solution hab oc exp(fci), and, therefore, unstable, of Einstein's equations around 
a super-extremal Kerr space-time. However, the explicit expression of the metric perturbation is so lengthy that it 
is hard to use it in further computations to evaluate, e.g., the perturbation effects at the ring singularity. One may 
want to compute the perturbed values of the Riemann invariants. A basis of algebraic invariants for the Riemann 



tensor in vacuum is given by the complex scalar fields [14 1 



W, := iCafccdC"'^'' = 2*0*4 - 8*1*3 + 6(*2)' 

W2 -ia6,dC^'^e/C"=-^"'' = 6*4*0*2 -6(*2)'-6(*i)2*4-6(*3)'*0 + 12*2*1*3 (13) 

o 

where Catcd the Weyl tensor and Cabcd ■= {C abed + i* C abed) (note that the Kretschmann invariant is given by the 
real part of Wi.) Since the only nonzero Weyl scalar for the Kerr background is *2 p=l, 15], the linearization of ([13]) 
yields 

5Wi = 12*2 (5*2, 5W2 = -18(*2)^ (5*2. (14) 

However, 5*2 — for arbitrary perturbations (with the exception of stationary, axially symmetric perturbations) of 
the Kerr spacetime jisj . and thus algebraic invariants are not modified to first order. A similar situation was found 
for the unstable modes of the negative mass Schwarzschild spacetime for which the perturbation effect on the 
singularity was then analyzed by computing differential invariants of the Riemann tensor, an approach that is hard 
to implement in this case, in view of the above mentioned length and complexity of the explicit expressions for the 
perturbed metric components. 
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We therefore at the moment do not know the effect of the perturbations constructed at the ring singularity and are 
not as confident as in the Reissner-Nordstrom case that the perturbations constructed are "conservative enough" in 
their behavior at the singularity. 

In any case, we must stress that the unstable (numerical) solutions of the Teukolsky equation found here are 
different from the algebraically special modes suggested in reference ■ In particular, the Starobinsky constant does 
not vanish, and, therefore, they represent a new type of solutions of the Teukolsky equation, of which we have found 
a few examples, through some simplifying assumptions, such as taking to = 0, etc. It would clearly be interesting to 
see what happens if these restrictions are lifted. We are currently working on this problem. 

Summarizing, we have shown explicitly in analytic form that the Reissner-Nordstrom space-time is linearly unstable 
when > M^, M > 0, even in the case the perturbations are "conservative" in the sense that they are small at the 
singularity in an appropriate sense. We have also numerical evidence that the Kerr space-time is unstable for a > M 
at least for some values of a, M . Further work is needed to confirm that the instabilities occur for all the range of 
parameters in super-extremality. 
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